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1. Introduction 



Lattice QCD is an important tool for the study of finite temperature QCD. Typically one is forced 
to use larger lattice spacing (a) for finite temperature quantities than zero temperature quantities. 
This is mainly because the temperature regime of interest Tc = {N-^a)^^ {N^ = 4 - 8 ) is about (1 
fm)^\ and simulations at several gauge couplings have to be performed to locate the transition 
temperature. 

This makes improved staggered fermion formulations attractive for thermodynamics studies. Im- 
provements of rotational symmetry via "knight's term" (p4 action) lai^gely reduces the cut-off ef- 
fects in the pressure and energy density at high temperatures ||l|]. Also, it is well known that 
staggered fermions break the flavor symmetry of continuum QCD and this effect can be largely re- 
duced if fat links are used in the fermion action |^. In this contribution we present an exploratory 
study of QCD phase diagram with a fermion action which combines the improvements in [|l|] and 
[|](p4fat7 action). A more detailed description of the action and simulation parameters are given 
in section ^. The phase transition for 3 degenerate flavors of p4fat7 quark is described in section ^. 
Comparison between p4fat7 and p4 action used in [^] is described in section ^ 



2. Action and simulation parameters 

For this study we use the tree-level improved Symanzik gauge action and p4fat7 action. 

S{x)=^Sg{x)+Sf{x), 
"^^W = L f 1 - ^ReTr[[/^W[/v(x + A)t/;(x + v)[/JW]) (2.1) 

^ 1 - ^ReTr [U^{x)U^{x + (L)Uy{x + 2ii)Ul{x + + v)ul{x + v)Ul{x)\ - ^ReTr[At ^ v] 



The fermion part of the action (5/) is shown in Fig. [T[ The coefficient in the smeared link are 
chosen to cancel flavor symmetry breaking at order a^a^ [Q], 

Here we use the above values of the fat7 coefficients without tadpole improvement. The effect of 
tadpole improvement on quenched configurations is studied in [Q]. 

All the codes needed for dynamical gauge evolution (HMD R algorithm) is implemented for QC- 
DOC in Columbia Physics System (CPS) and numerically checked against existing Asqtad and p4 
action results. The "knight's move" term of p4 action is implemented without non-nearest commu- 
nication by breaking the term into Parallel transport ( ^'^{x) = U^{x)y{x + fi)) and recombination 
[R^ (x) = Lv^^i ( V^(-^ + v) ± y{x - v)) , as illustrated in Fig. g 

Inversion of p4fat7 Dirac operator, which is the dominating part of gauge evolution for small 
masses, is currently running at '-^ 31% of the peak for 4^ x 6 local volume, which is the local 
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Figure 1: Diagram of p4fat7 action. The numbers above graphs denote the muhipUcities of each graphs. 
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Figure 2: Illustration of implementation of p4 action. The arrows and circles denote the parallel transports 
and recombinations respectively. 



volume for 32^ x 6 on a 5 12-node QCDOC partition. Further optimization of the Dirac operator, 
especially at smaller volumes, is possible. Rational Hybrid Monte Carlo(RHMC)[^ implementa- 
tion is also in progress. All the simulations were done with R algorithm[|6|] with the step size equal 
to 40% of the bare mass in lattice units. 



3. Nf = 3 phase transition with p4fat7 

Fig. m shows the coupling at the transition as a function of quark masses in lattice units. j8c are 
located by measuring susceptibilities in average plaquette, Polyakov loop and chiral condensate. 
The location of peaks coincide within errors. The use of tree-level coefficients enables the use of 
the reweighting procedure, which is used in locating the transition temperature. At /3c for each 
mass, separate measurements of heavy quark potential and spectrum calculations were done on 
16^ X 32 lattices for physical scale setting. 
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Volume 


sites/node 


CG(MFlops/node) 


Smearing(MFlops/node) 


2242 


64 


173 


271 


44 


256 


207 


313 


4^x62 


576 


264 


264 


64 


1296 


231 


164 



Table 1: Performance of p4fat7 evolution codes on QCDOC, 420Mhz, 1024 nodes, MFlops/node 
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Figure 3: Transition coupling in lattice units 
(j3e) for 3 flavor p4fat7 action. The curves rep- 
resent fits to j3c(ffJ^) = A -\- B{mga)'^ . Red and 
green points denote Nr = 4 and Nr ~ 6 respec- 
tively. 




Figure 4: The susceptibility in chiral condensate 
for 8^ X 4 



It has been shown that staggered type fermions exhibits a transition from a rapid crossover to a 1st 
order transition in small masses. Initial studies suggest that the mass where the transition occurs 
for Nt: = 4 may be larger than [^], as shown by the rapidly increasing peak in susceptibility in 
Fig. ^ This may indicate the additional smearing does not suppress finite lattice spacing error 
much. Further study is in progress. 

Fig. ^ shows the heavy quark potential and the lattice spacing, measured by ro|^. The error in each 
r^ro includes the difference between 2 different spacial smearings. 



4. Comparison with p4fat3 action 



P4fat3 action, used in |^, has the knight's move term, 1-link term as well as 3-staple terms. In the 
notation in Fig. |l|, the coefficients are 

3 1 3 CO 

^^ = 4TT6^' ^^ = 4TT6^' ^^=^^ = 0' « = 0-2- (4-1) 
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Figure 5: Scale from heavy quark potential for 8x4 p4fat7 lattices 
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Figure 6: Comparison of Tc/my between p4fat7 and the p4fat3. 



Fig. g shows the behavior of the transition temperature for the p4fat7 and p4fat3 actions. The 
transition temperature in units of vector meson mass shows a good agreement between the 2 actions. 
A further study will be needed for the continuum extrapolation. 



5. Summary and Future plans 

We studied the phase transition of p4fat7 action for 3 degenerate flavors of quarks. p4fat7 action 
exhibits the finite temperature phase transition and appears to be suitable for other studies such as 
Equation of States. Although the use of fat7 smearing shows the expected reduction of flavor sym- 
metry breaking symmetry breaking in smaller lattice spacings[Q], we do not see a large improve- 
ment in flavor symmetry breaking for the lattice spacing relevant for A^^- = 4. More improvement 
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such as employing tadpole improvement may be needed. We also showed some preliminary results 
for Nx = 6 phase transition. More systematic study of Nf = 3 and 2+1 transition is in progress. 
Dirac operator and other routines needed for p4fat7 Hybrid Moleculardynamics are optimized for 
QCDOC and performing at ~ 30% for the volumes relevant for finite temperature studies. Imple- 
mentation for RHMC is under way. 
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